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I. INTRODUCTION 

Recent experiments have discovered^ that the resis- 
tivity of a high-mobihty two-dimensional electron gas 
(2DEG) in GaAs/AlGaAs heterostructures subjected to 
microwave radiation of frequency w exhibits magne- 
tooscillations governed by the ratio uj/oJc, where tOc is 
the cyclotron frequency. Subsequent work^~^ has shown 
that for samples with a very high mobility and for high 
radiation power the minima of the oscillations evolve into 
zero-resistance states (ZRS). 

These spectacular observations have attracted much 
theoretical interest. As was shown in Ref. 8, the ZRS can 
be understood as a direct consequence of the oscillatory 
photoconductivity (OPC), provided that the latter may 
become negative. A negative value of the OPC signifies 
an instability leading to the formation of spontaneous- 
current domains showing zero value of the observable re- 
sistance. Therefore, the identification of the microscopic 
mechanism of the OPC appears to be the key question 
in the interpretation of the data.^""^ 

A mechanism of the OPC proposed in Ref. 9 is based 
on the effect of microwave radiation on electron scat- 
tering by impurities in a strong magnetic field (see also 
Ref. 10 for an earlier theory and Ref. 11 for a systematic 
theory). An alternative mechanism of the OPC was re- 
cently proposed in Ref. 12. In contrast to Refs. 9-11, this 
mechanism is governed by a radiation-induced change of 
the electron distribution function. Because of the oscilla- 
tions of the density of states (DOS), v{£), related to the 
Landau quantization, the correction to the distribution 
function acquires an oscillatory structure as well. This 
generates a contribution to the dc conductivity which os- 
cillates with varying LojiOc- A distinctive feature of the 
contribution of Ref. 12 is that it is proportional to the 
inelastic relaxation time Tin. A comparison of the re- 
sults of Refs. 11 and 12 shows that the latter contribu- 
tion dominates if Tin Tq (where Tq is the quantum, 
or single-particle, relaxation time due to impurity scat- 
tering), which is the case for the experimentally relevant 



temperatures. 

The consideration of Ref. 12 is restricted to the regime 
which is linear in both the ac power and the dc electric 
field. The purpose of this paper is to develop a com- 
plete theory of the OPC governed by this mechanism, in- 
cluding nonlinear effects. We will demonstrate that the 
conductivity at a minimum becomes negative for a large 
microwave power and that a positive sign is restored for 
a strong dc bias, as it was assumed in Ref. 8. 

The paper is organized as follows. First, in Sec. II we 
formulate a general approach to the problem. In Sec. Ill 
we calculate the non-equilibrium distribution function for 
overlapping Landau levels (LLs). In Sec. IV we consider 
the OPC in the linear regime with respect to the dc field. 
In Sec. V we analyze the ZRS and calculate the sponta- 
neous electric field in the domains. In Sec. VI we turn 
to separated LLs. Section VII deals with the inelastic 
relaxation due to electron-electron scattering. Finally, 
in Sec. VIII we briefly discuss the magnetooscillations in 
the Hall photoresistivity. In Sec. IX we summarize our 
results and compare them with the experimental data. A 
brief account of the results of this paper was presented 
in Ref. 13. 



II. GENERAL FORMALISM 

We consider a 2DEG (mass to, density rig, Fermi ve- 
locity vf) subjected to a transverse magnetic field B = 
{rnc/e) lOc- We assume that the field is classically strong, 
^cTtt 2> 1, where Ttr is the transport relaxation time at 
B The photoconductivity Cph determines the longi- 
tudinal current flowing in response to a dc electric field 
^dc, j ■ fdc = fph^dc' ™^ presence of a microwave elec- 
tric field £uj cosLot. The more frequently measured"'^"'*'^''' 
longitudinal resistivity, pph, is given by pph ~ plyC^ph, 
where pxy — eB /n^c is the Hall resistivity, affected only 
weakly by the radiation. 

We start with the formula for the dc conductivity per 
spin: 
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0-ph = j deadcie) [-d^fie)], 



(1) 



where /(e) is the electron distribution function, and 
adc{s) determines the contribution of electrons with en- 
ergy £ to the dissipative transport. In the leading 
approximation/^'^^ 



o-dc(£) 



^tr,B(£) 



tr,B 



(2) 



where Tir B is the transport scattering time in a quantiz- 
ing magnetic field, Ttr,B(£) = Ttiiyo/i'{e), and uq = m/27r 
is the DOS per spin at zero B (we use h = 1). We note 
that Eq. (2) has a Drude form with the DOS i^{e) and 
the transport time Ttr,B(£) dependent of energy due to 
the Landau quantization. 

For a classically strong magnetic field, cOcTu ^ 1, the 
above expression reduces to 



(Jdcie) = C7dci> (e) 



(3) 



where = e^t'oUp/2ti)^rtr is the dc Drude conductivity 
per spin in strong B and we introduced the dimensionless 
DOS, i>{s) = u{e)lua. 

We neglect here the effect of the microwaves on the im- 
purity collision integral, which yields a subleading con- 
tribution to the pho conductivity, as discussed in Sec. IV. 
The dominant effect is due to a non-trivial energy depen- 
dence of the non-equilibrium distribution function /(e). 
The latter is found as a solution of the stationary kinetic 
equation for the zero angular harmonic of the distribu- 
tion function /(e): 



St^{/} + Stdc{/} = -Sti„{/} , 



(4) 



Here the left-hand side represents the effect of the mi- 
crowaves (St„) and of the dc field (Stdc) in the presence 
of impurities while the right-hand side accounts for the 
inelastic relaxation. 

The first term on the left-hand side describes the ab- 
sorption and emission of microwave quanta; the rate of 
these transitions was calculated in Ref. 12, yielding 



St,{/} 



^tr,B 



i(e + u;)[/(e + u;)-/(e): 



4a;2 ^ 2(u; ± uJc? + T-%{e + w) + T-%{e) 



+ {w -w} . 



(5) 



We will assume that |tj ± ^ '^tTs' ^^^^ excluding a 
narrow vicinity of the cyclotron resonance. This allows 
us to neglect the e dependent terms in the denominator, 
which reduces Eq. (5) to the form 



Sta.{/} = Si 



E 

± 



where the ac Drude conductivity per spin is given by 



4Ttr(w±Wc)^' 



(7) 



Furthermore, we assume, in accordance with the exper- 
iments, a linear polarization of the microwaves. For a 
circular polarization, one should retain only one term 
on the right-hand side of Eq. (7), which, away from the 
cyclotron resonance, does not affect the results in any 
essential way. 

The second term on the left-hand side of Eq. (4) rep- 
resents the effect of the dc electric field. The impu- 
rity scattering in a quantizing magnetic field leads to 
the spatial diffusion with a diffusion coefficient D^{e) = 
t;|./2a;^rtr,B(£)- In view of conservation of the total en- 
ergy e + eSdcX in a dc field (directed along the x axis) 
the spatial diffusion is translated into the diffusion in e- 
space, with a diffusion coefficient {eSdcY' D-Q{e) and the 
DOS v{e), 



Stdc{/} = y-\e) de [v{e) e'^SiMs) /(e) ; 



(8) 



Equation (8) can also be obtained from Eq. (6) by taking 
the limit a; ^ and replacing the period-average of the 



ac field squared, 2^,^, by £^^. 

Substituting Eqs. (6), (8) in the kinetic equation (4), 
we get 



Sl^^Y.^{e±u;)[f{e±u;)- f{e)] 



± 



d,[u^{e)d,f{e)] = 



m-Me) 



(9) 



where the inelastic processes are included in the relax- 
ation time approximation and /t(£) is the Fermi distri- 
bution. A detailed discussion of the inelastic relaxation 
and a calculation of the inelastic relaxation time Tin are 
relegated to Sec. VII. 

Equation (9) suggests convenient dimensionless units 
for the strength of the ac and dc fields: 



Tin 
Ttr 



Qdc = 



2 Tin / eSdcVp 
Ttr \ UJc 



(a;2_^2)2 

TT 



(10a) 
(10b) 



With these notations, Eq. (9) reads 
^5^i>(e±a;)[/(e±a;)-/(e)] 



%±w)[/(e±u;)-/(e)], (6) ^4^2^(e) 



d,[u\6)dJis)]=fis)-Me). (11) 



Note that and Qdc are proportional to Tin and are in- 
finite in the absence of the inelastic relaxation processes. 
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The left-hand side of the kinetic equation (11), as well 
as Eqs. (1),(3) for the photoconductivity, can also be ex- 
tracted from the quantum Boltzmann equation of Ref . 1 1 , 
as we show in Appendix. Below, we use Eqs. (1),(3),(11) 
to analyze the photoconductivity in the both limiting 
cases of overlapping and separated LLs. 



III. NON-EQUILIBRIUM DISTRIBUTION 
FUNCTION INDUCED BY MICROWAVE 
RADIATION 



We consider first the case of overlapping LLs, with the 
DOS given by 

27re / TT \ 
i>=l-2(5cos , (5 = exp < 1. (12) 

Here Tq is the zevo-B single-particle relaxation time, 
which is much shorter than the transport time in high- 
mobility structures, Tq <C rtr (because of the smooth 
character of a random potential of remote donors). The 
existence of a small parameter S simplifies solution of the 
kinetic equation (11). To first order in S, we look for a 
solution in the form 



I I I 




(e - ep)/cOj, 



FIG. 1. Schematic behavior of the oscillatory density of 
states !/(e) and radiation induced oscillations in the distribu- 
tion function /(e) for sm{2nLj/ujc) > 0. 

One might naively think that the small correction (14) 
to the distribution function will only weakly affect the 
conductivity. This is not the case, however. The reason 
is that, due to the fast oscillations in /osc, the deriva- 
tive 9e/osc may be large. As a result, a small variation of 
the distribution function (14) can strongly affect the con- 
ductivity, Eq. (1). In particular, when the regions of an 
inverted population in /(e) correspond to the maxima in 
i'(e) (as in Fig. 1), the linear-response conductivity may 
become negative, as we show below. 

IV. OSCILLATORY PHOTOCONDUCTIVITY 



/ = /O + /osc + 0(<5'), fo: 



6 Re 



(13) 



We assume that the temperature (measured in energy 
units, fee = 1) is sufficiently high, T 3> ujc, implying a 
scale separation between the smooth energy dependence 
of functions fo,i{s) on a scale of the order of T and the 
fast oscillations with a period ojc We also assume that 
the electric fields are not too strong [^'(^(w/T)^ ^ 1 and 
Qdc(wc/T)^ <C 1], so that the smooth part /0(e) is close 
to the Fermi distribution fris) at a bath temperature T; 
otherwise, the temperature of the electron gas is further 
increased due to heating. We obtain^^ 



, , , c^c 9/t . 27re sin ^ + 4Qd, 

— ^ 7: TT- sm 



sin2 ^ + Qdc 



(14) 



Thus, the oscillations of the DOS v{e) induce an oscilla- 
tory contribution /osc(£) to the distribution function, as 
illustrated in Fig. 1. 



To calculate the photoconductivity, we substitute 
Eq. (14) for the distribution function into Eq. (1). Per- 
forming the energy integration in Eq. (1), we assume (in 
conformity with the experiment) that T is much larger 
than the Dingle temperature, T ^ l/27rTg. Under this 
condition, the temperature smearing yields a dominant 
damping factor of the Shubnikov-de Haas oscillations, 
X/sinhX with X = ^tP-TIu:^- The terms of order 8 
in Eq. (1) are then exponentially suppressed, 

d de cos OJt oc (5exp(-27r^T/wc) < (5^, 

J i^c 

and can be neglected. The leading w-dependent contri- 
bution to (Tph comes from the term generated by the 
product of c?e/osc(£) (5cos^ and the oscillatory part 
—25 cos of £'(e). This term does survive the energy 
averaging, — J de cos^ d^fx — 1/2. We thus find 



1 + 25^ 



27TUJ , 



4Qdc 



1 



Qdc 



(15) 



Equation (15) is our central result. It describes the 
photoconductivity in the regime of overlapping LLs, in- 
cluding all non-linear (in S,^ and £dc) effects. Let us 
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analyze it in more detail. In the linear-response regime 
(^dc ~^ 0) and for a not too strong microwave field, 
Eq. (15) yields a correction to the dark dc conductiv- 
ity (Tdc = <^dci^ + 25^), which is linear in the microwave 
power: 



gph — gdc 



TTU) 27ru> 
— sin 



(16) 



in agreement with Ref. 12. It is enlightening to com- 
pare Eq. (16) with the contribution of the effect of the ac 
field on the impurity scattering. The analytic result, 
Eq. (6.11) of Ref. 11, in the notation of Eq. (10) is 



Jii] 



fdc 



CTdc 



-12-3-^2^^ — sin -h sin^ — 



(17) 



This result has a similar frequency dependence as 
Eq. (16); however, its amplitude is much smaller at 
Tin ^ Tq, i.e., the mechanism of Refs. 9-11 appears 
to be irrelevant. Physically, the effect of the ac field 
on the distribution function is dominant because it is 
accumulated during a diffusive process of duration Tin, 
whereas Refs. 9-11 consider only one scattering event. 
Apart from the magnitude, the two contributions are 
qualitatively different in their temperature and polariza- 
tion dependence. Specifically, the contribution related 
to the change of the distribution function is strongly 
temperature-dependent (due to the T-dependence of Tin, 
see Sec. VII) and does not depend on the direction of the 
linear polarization of the microwave field. On the other 
hand, the effect of microwaves on the impurity collision 
integral yields a T-independent contribution which de- 
pends essentially on the relative orientation of the fields 
£^ and fdc [Eq. (17) represents the result averaged over 
the polarization direction] . 

With increasing microwave power, the photoconduc- 
tivity saturates at the value 



^ = l_852![^^cot^ 



Vu^sir? — > 1. 



(18) 



Note that although the correction is proportional to 
5'^ <^ 1, the factor 8tt{uj /cUc) cot{TTUJ /ujc) is large in the 
vicinity of the cyclotron resonance harmonics lu = ktOc 
{k = 1,2,...), and allows the photo- induced correc- 
tion to exceed in magnitude the dark conductivity cTdc- 
In particular, dph around minima becomes negative at 
> P^ > 0) with the threshold value given according 
to Eq. (15) by 



V,., = I Ad — sm sm — 



(19) 



The evolution of the B dependence of the photoresistiv- 

ity with increasing microwave power is illustrated in 
Fig. 2. More specifically, the curves in Fig. 2 correspond 
to different values of the dimensionless parameter 



pW=P^(a;e = 0) 



e Tin SttP 



Tie Ttr W^S hu}"^ ' 



(20) 



where P = c\£u\^S/8'k is the microwave power over the 
sample area S, c is the speed of light, and we restored 
the Planck constant for convenience. 
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FIG. 2. Photoresistivity (normalized to the dark Drude 
value) for overlapping Landau levels vs uJc/oj at fixed 
LJTq = 2n. The curves correspond to different levels of mi- 
crowave power P^ ^ — {0.24, 0.8, 2.4}. Nonlinear I-V char- 
acteristics at the marked minima are shown in Fig. 3. 



V. ZERO-RESISTANCE STATES 

Let us now fix lo/lOc such that P* > 0, and consider 
the dependence of Uph on the dc field S^c at a suffi- 
ciently strong microwave power P^^ > P* , corresponding 
to the negative linear-response photoconductivity. As 
follows from Eq. (15), in the limit of large £dc the con- 
ductivity is close to the Drude value and thus positive, 
fph = (1 ^ 6(5^)(T^^ > 0. Therefore, cTph changes sign at 
a certain value of the dc field (sec Fig. 3), which is 
determined by the condition Qdc = {Plu - P^)/P^- The 
negative-conductivity state at S^c < unstable with 

respect to the formation of domains with a spontaneous 
electric field of the magnitude £^^.^ Using Eqs. (10), (19), 
we obtain 



£* = ^ — ( — 



1/2 



= - 



- 1 

-,1/2 



1/2 



2a;2(w2 



X -Re 46 — sin sm — 



1/2 



(21) 



with £* being the threshold value of the ac field at which 

the ZRS develops and = vp / lUc the cyclotron radius. 
Equation (21) relates the electric field formed in the do- 
mains (measurable by local probe^) with the excess power 
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of microwave radiation. It is worth noticing that the last 
expression for does not exphcitly contain the rate of 
the inelastic processes. 




FIG. 3. Current-voltage characteristics [ dimensionless 
current jx = ((7ph/crJ;.)£dc vs dimensionless field £dc = Q^i^] 
at the points of minima marked by the circles in Fig. 1. The 
arrows show the dc field in spontaneously formed domains. 



VI. SEPARATED LANDAU LEVELS 



We now turn to the regime of strong B, uicTq/n » 1, 
where the LLs get separated. The DOS is then given 
(within the self-consistent Born approximation) by a se- 
quence of semicircles of width 2r = 2(2ct;c/7''Tq)^/^: 



r2 - (e - 



■0Jc/2f. (22) 



We use Eqs. (1) and (11) to evaluate the OPC at Qdc 
to first order in and estimate the correction of the sec- 
ond order. The condition T ^ Uc allows us to separate 
the slow dependence on e on the scale of T and fast oscil- 
lations with the period Uc in the integral (1) by averaging 
over the period of the oscillations. After integrating the 
resulting slow-varying functions we obtain 



'dc 



{[v{e + uj) - v{e ^ oj)]d,v\e)) 



4 

+ '^^w Xlfe I ak,i{v{£ + kLo)i'{£ + li^) d^iP- (e) )^ 



(23) 

Here the angular brackets denote averaging over e within 
the period Wc, and a^, ; are numerical coefficients. The 
result for separated LLs, Eq. (22), reads 

16wc 



CTph 



STT^r 



$(a;) = — Re 
47r 



1-Pc 



o 



(24) 



arccos(|a;| — 1) — 



VN(2-N) 



The photoresistivity for the case of separated LLs, 
Eq. (24), is shown in Fig. 4 for several values V^^ of the 
microwave power. Notice that a correction to Eq. (24) of 
second order in is still small even at microwave power 
exceeding the threshold value 



(25) 



since ujcVl^jV ^ r/o; <C 1. This means that it suffices to 
keep the linear-in-T^^^ term only even for the microwave 
power Vt^ > V* at which the linear-response resistance 
becomes negative. 




CO /CO 



FIG. 4. Photoresistivity (normalized to the Drude value) 
for separated Landau levels vs ujc/w at fixed WTq = 16tt. 
The curves correspond to different levels of microwave power 
VL°^ = {0.01, 0.03, 0.05}. 

As in the case of overlapping LLs, a negative value 
of the linear-response conductivity signals an instability 
leading to the formation of domains with the field S^c 
at which crph(fdc) ==0. It turns out, however, that for 
separated LLs the kinetic equation in the form of Eq. (9) 
yields zero (rather than expected positive) conductivity 
in the limit of strong £dc-^^ This happens because elastic 
impurity scattering between LLs, inclined in a strong dc 
field; is not included in Eq. (9). The intcr-LL transitions 
become efficient in dc fields as strong as [see Eq. (5.5) of 
Ref. 11] 



dc 



Ttr 



1/2 



(26) 



which actually gives the strength of the field in domains 
for the case of separated LLs. 



VII. INELASTIC RELAXATION DUE TO 
ELECTRON-ELECTRON COLLISIONS 



Finally, we calculate the inelastic relaxation time Tin- 
Of particular importance is its T dependence which in 
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turn determines that of CTph. At not too high T, the dom- 
inant mechanism of inelastic scattering is due to electron- 
electron (e-e) collisions. It is worth emphasizing that the 
e-e scattering does not yield relaxation of the total en- 
ergy of the 2DEG and as such cannot establish a steady- 
state dc photoconductivity. That is to say the smearing 
of fo[e) in Eq. (13), which is a measure of the degree of 
heating, is governed by electron-phonon scattering. How- 
ever, the e-e scattering at T 3> t^^c does lead to relaxation 
of the oscillatory term /osc, Eq. (14), and thus determines 
the T behavior of the oscillatory contribution to CTph- 

Quantitatively, the effect of electron-electron interac- 
tion is taken into account by replacing the right-hand 
side of Eq. (9) by — Stee {/}, where the collision integral 
Stee {/} is given by 

Stee{/} = j de' j dEA{E, e, e') 

-/(£)/W(£+)/(e')/^'H£-) 

+/('^)(e)/(£+)/('')(e')/(e'-)] , (27) 

and f'He) = 1 - fie), e+ = e + E, e'_ = e' - E. The 
function A{E, e, e') describes the dependence of the ma- 
trix element of the screened Coulomb interaction on the 
transferred energy E and the energies of colliding parti- 
cles. In what follows we use the linearized form of Stee- 



A. Overlapping Landau levels 

For overlapping LLs, we put i/ = 1 in accord with the 
accuracy of Eq. (15). Then the kernel in Eq. (27) de- 
pends on the transferred energy E only and is given by 
the general formula^^ 



d^q 



|C/|2(Re(P))2 



(28) 



Here the factor of two accounts for spin, {T>){E, q) is the 
angle-averaged particle-hole propagator, 



U{E,q) = 



q + K{l + iE{V)) 



(29) 



the dynamically screened Coulomb potential, and k = 

Ane^h'Q the inverse screening length. 

The propagator {T>){E,q) has the (disorder- 
independent) Fermi-liquid form for ultra ballistic mo- 
menta, q^ qi = (wc/D)^/^, where D = i?^/2Ttr is the 
diffusion coefficient in a classically strong magnetic field. 



{V) = {q'v'p-E') 



-1/2 



g » gi = {uJc/D) 



1/2 



(30) 



With lowering q it crosses over into the quasiclassical 
particle-hole propagator in a smooth random potential 
("ballistic diffuson")i9'20 



n=0 



-i{E-nuJc) + Dq^ + in?lTtr' 



q^qi 



(31) 



where Jn{x) is the Bessel function. 

Let us calculate the kernel A{E) for E within the 
(experimentally relevant) temperature range, 



(32) 



The integration domain in Eq. (28) naturally divides into 
two parts: q > qi and q < qi. In view of Eq.(32) the 
propagator (30) for ultra-ballistic momenta q^ qi does 
not depend on the transferred energy E ^ T; specif- 
ically, (V) ~ l/qvp, while the screening is effectively 
static, U{E,q) ~ U{0,q). It follows that the contribu- 
tion A-^{E) oi q> qi to the integral (28) has the form 



A>iE) = 



2-^5 F qi 



(33) 



By contrast, the contribution of g < to Eq. (28) is 
strongly _E-dcpendcnt and is found as a sum over peaks 
coming from different n in Eq. (31), 



A<iE) = Y,^:;{E)- 



(34) 



When calculating the contribution of a single peak 
A^, the following approximations are justified: in 
Eq. (31), the product of the Bessel functions J^{qRc) 
for the relevant momenta q » can be replaced by 
2cos^ ipn/nqRc with ipn = qRc + 7r(2n + l)/4; the term 
n^/rtr in Eq. (31) can be neglected at E ^ T in view 
of Eq. (32); also, the first term in the denominator of 
Eq. (29) can be omitted for g < gi <C k. Introducing the 
dimensionless parameters An{E) = \E/ujc — n\ < 1/2, 
Pn = 27r2wcTtr/n2 » 1, and jniE) = pl'^AniE) we ex- 
press the contribution of the nth peak for A„ <IC 1 as^^ 



0' 



f 

Jo 



dx 



4cos^(a„a;) 



+ 1 



(35) 



After averaging over the fast oscillations with UnX, where 

1 /3 

an = nPn /ti", the integration yields 



1 [ c0\ A„<c/3-^/' .gg. 



where c = r(7/6)/2i/3(37r)i/2r(2/3) ~ 0.18. We now use 
Eqs. (27), (33)-(36) to calculate the relaxation rate t^~^ 
of the oscillatory component of the distribution function, 
./osc(e)- As we are going to show, it is dominated by the 
large-momentum transfers, q ^ qi. The contribution 
of this large-momentum region is easily evaluated: since 
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A^{E) is energy- independent, it is sufficient to tafce into 
account the out-scattering term only. We thus return to 
the right-hand side of Eq. (9) with the inelastic relaxation 
rate 7^~^ replaced by^^ 

oo 

= [dE A>{E) f ( coth ^ - tanh 



■In 



A-KEF Wc(WcTtr)^/^ 



(37) 



Note that in Eq. (37) depends on B through the 
logarithmic factor only and crosses over into the conven- 
tional zero-B result t^^^ = {nT'^/AeF) ln{KVF/T) when T 

exceeds aJc('^'cTtr)"'"^^- 

Let us now turn to the contribution of the region 
q ^ qi- Evaluating the out-scattering term [similar to 

the first line of Eq. (37), with A>{E) replaced by A<{E), 
Eqs. (34), (36)], we find (t<)-i ~ {cOcTtry^^ /ef, which 
can exceed the large-g contribution (37). However, the 
relaxation time approximation is no longer valid for q < 
qi. Indeed, the main contribution to {t^}^^ comes from 

the energy transfers E close to nuc, A„(£') 0n^^^ , 
see Eq. (34). Such processes are inefficient as far as the 
relaxation of /osc is concerned, since the energy transfer 
is almost commensurate with the period of the oscilla- 
tions in the distribution function. In other words, the 
out scattering term is almost compensated by the in- 
scattering one, so that A^{E) = A'!^{E) is effectively 
replaced by '^^^ A'^A^{E). As a result, the contribution 
of region of q < qi to the relaxation rate is dominated by 
q ^ qi and is given by Eq. (37) without the logarithmic 
factor, and thus can be neglected. 

In fact, the situation is similar to the momentum relax- 
ation due to small-angle scattering off a smooth random 
potential. The momentum relaxation time Ttr and the 
out-scattering (or single-particle) relaxation time Tq in 
that case differ by the factor (1 — cost/)), which accounts 
for a reduced contribution to the resistivity of small-angle 
scattering with <C 1: 



^Aw{2kFSu.t)^[] 

2-K ^ 2' 1 (1 - cos^) 



where W{q) is the Fourier transform of the correlation 
function of the random potential. A similar result for 
the relaxation of the oscillatory part of the distribution 
functions is obtained from Eq. (27). We linearize Eq. (27) 
with the distribution function in the form (we assume 
T>a;) 



/ = /t + vie) d,h 



(38) 



where ^p{e) = ip{e + Wc)- This ansatz is suggested by 
Eq. (14) and will be confirmed by the calculation below. 
We substitute Eq. (38) in Eq. (27) and use the condition 
T ^ (jjc which allows us to separate the slow dependence 
on E, s' on the scale of T and fast oscillations with the 



period Uc by averaging over the period of the oscillations. 
Using A{E) = A{E + uJc) [see Eqs. (33), (34), and (36)], 
we obtain the following integrals over the slow variables 
which all produce the same result, 

-Jde'J dE[f^'\e+)Ms')f^'\e'_) 

+ fT{e+)f^^\e')fT{e'_)WT{e) 
= jde' j dE [/T(e)/T(£')4'^(e-) 

+ f!^\e)f!p\e')fT{e'_)]dMe+) 
= -Jds'JdE [fT{e)f^^\e+)fP{e'_) 

+ f^^\e)fT{e+)fT{e'M'fT{e') 
= jde' j dE [fT{e)f!^\e+)fT{£') 



deMe) . 



(39) 



The collision integral thus reads 

-Stin {/} = 



2 Be 

X ( A{E) [ ^(£) - ^(£ + E)+ ip{e') -^{e' -E)]),,,E , 

where the angular brackets denote averaging over e' and 
E within the period ujc- For a harmonic modulation of 
the distribution function, ip{e) oc cos(27r£/a;c + 0), as in 
Eq. (13), and using A{E) = A{-E), we obtain 

-Stin{/}=^^^, 
Tin 



-{A{E)[l-cos{27rE/uJc)])E. (41) 

Because of the factor 1 — cos(27r£'/u;c) the contribution of 
small momenta (34), (36) to the relaxation rate is small 
compared to that of g > gi. In the latter case, due to 
the energy- independent kernel A(E), Eq. (33), the in- 
scattering part is zero on average, and 7^~^ coincides with 
the out-scattering rate (37). 

The inelastic relaxation time Tin as obtained above 
[Eqs. (37), (41)] depends on energy e.^^ This makes the 
problem somewhat more complicated than the model 
considered in Sec. II with a phenomenological, e- 
independent parameter Tin- However, characteristic ener- 
gies are e ^ T [see Eq. (38)], so that the £-dependence in 
Eqs. (37), (41) does not change the T"^ scaling of Tin but 
only yields a numerical factor. In particular, repeating 
the analysis of Eq. (9) in the linear in Vuj regime with 
the £-dependent Tin, we find that Tin entering Eq. (10a) 
is effectively replaced by 

Jden4e,T){-defT) 

= Tin(0, T) I rf£^-Z^|g_ ^ 0.822 Tin(0, T). (42) 



7 



Therefore, at T a; the linear photoresistivity, Eq. (16), 
scales as T^^. 

At T <i; w (but still T » Wc), the e-dependence in Tin 
becomes more important. Indeed, solving Eq. (9) to the 
linear order in E"^ (and at fdc = 0), we find the following 
oscillatory contribution to the distribution function 

X ^%±a;)[/T(e±a;)-/T(e)]. (43) 

± 

Equation (1) then reproduces Eq. (16) for the photocon- 
ductivity, with Tin (entering Vuj) given by 

Tin = j ^ri„(£, T)[fT{e - a;) - /t(£ + a;)]. (44) 

For T ^ CO this expression reduces back to Eq. (42), while 
in the opposite limit T <C w it yields 

Tin = I ^Ti„(£,T) = ^nn{0,T). (45) 

Thus, the scaling of the photoresistivity at T a; 
transforms into the behavior for T <C w. 



B. Separated Landau levels 

We turn now to the case of separated LLs. Let us 
assume again that T is not too high, 

We « T < r(Ttr/Tq)l/' . (46) 

The dominant contribution, similarly to the case of 
overlapping LLs, is given by large transferred momenta 
q > (^/^)B)^/^ where Db = Rl/2Tt,^B is the diffusion 
coefficient in a quantizing magnetic field (see Sec. II). 
However, the situation for a strongly oscillating DOS, 
Eq. (22), is different in that it is no longer sufficient 
to deal solely with the out-scattering processes even for 
large g; that is the whole collision integral (27) should 
be taken into account. The kernel in Eq. (27) may be 
re-written as: 

A{E, e, e') 

where the function 11 for large q {T/DbY^^ reads^^ 
ILl^, = TTVQ{qvF)-^v{e)i>{e') . (48) 

The procedure leading to Eq. (40) is applicable in the 

case of separated LLs as well. As a result, the inelastic 
collision term on the right-hand side of Eq. (9) should be 
replaced by 



X {i){e)i}{e + E)i>{e')i){e' - E) 

X [ ifis) - ip{s + E) + ip{e') - ^{e' - E) ] ),,,e ■ (49) 

In other words, Eq. (9) becomes an integral equation for 
the periodic function (p{s) characterizing oscillations of 
the distribution function, 

Al^iE + Lj) - i>{e - oj)] = { u{s + E) y{e') y{e' - E) 
X [ip{e) - ip{e + E) + ^{e') - ^{e' - E)\ ),,,e , (50) 

where ^ is a smooth function of e, 

Analytical solution of Eq. (50) does not seem feasible. 
However, up to a factor of order unity,^^ we can rewrite 
the exact collision integral (49) in the relaxation-time ap- 
proximation, thus returning to Eq. (9) with 

One sees that the T and e dependence of Tin in the 
regime of separated LLs is the same as for overlapping 
LLs [Eqs. (37), (41)]. Therefore the temperature scaling 
of the linear- in- ■Pt^ photoresistivity, Eq. (24), is the same 
as found in Sec. VII A. In particular, cjph — Cdc scales as 
T-2 for T > w. 



VIII. OSCILLATORY HALL RESISTIVITY 

Finally, let us briefly discuss the issue of the microwave 
induced cj/wc^oscillations of the Hall resistivity p^y (an- 
tisymmetric part of pxy, i.e. = —Py^) detected in re- 
cent experiments.^^ The experimentally observed oscilla- 
tions of Pxy demonstrate the following properties. First, 
they have the same period and the opposite phase as 
compared to the oscillations of the dissipative resistivity 
Pxx- Second, the amplitude of the oscillations in pxy is 
roughly the same as in pxx- Third, the radiation-induced 
contribution 6pxy is odd with respect to the magnetic 
field B, which implies that the measured dpxy is indeed 
a contribution to the Hall resistivity pxy These observa- 
tions cannot be explained within either the mechanism 
related to the effect of the microwaves on the distribu- 
tion function (Ref. 12 and this work) or the one related to 
the effect on the elastic collision integral (Refs. 9-11), if 
one assumes, as usual, that the electron density rig is con- 
stant. If Ue was constant, the leading odd-in-B correction 
S(Ty^ to the Hall conductivity tr^, = ceUe/B should be 
smaller than 6axx by a factor l/wcTtr, which is of order 
10~^ under the experimental conditions. Therefore, the 
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observed Sp^y^ although very smaU in comparison with 
the bare pxy, appears to be two order of magnitude larger 
than what one would expect from the theory of the os- 
cillatory contribution to Pxy neglecting the oscillations of 
rie.^^ A possible origin of the observed oscillatory p^y may 
be a weak variation of Ue with lojujc- The observation^* 
of a variation of the period of the Shubnikov-de Haas os- 
cillations, ^Uxx oc cos{cnc/eB), appears to support the 
idea of the microwave-induced oscillations of the electron 
density. The issue warrants further study. 



IX. CONCLUSION 

To summarize, we have presented a theory of magne- 
tooscillations in the photoconductivity of a 2DEG. The 
parametrically largest contribution to the effect is gov- 
erned by the microwave- induced change in the distribu- 
tion function. We have analyzed the nonlinearity with 
respect to both the microwave and dc fields. The result 
takes an especially simple form in the regime of over- 
lapping LLs, Eq. (15). We have shown that the magni- 
tude of the effect is governed by the inelastic relaxation 
time (44), (41), (51), and increases as or (de- 
pending on the relation between T and to) with lowering 
temperature. For a sufficiently strong microwave power 
the linear-in-dc-field photoconductivity becomes negative 
leading to formation of domains with zero resistivity. We 
have calculated the threshold power at which this zero- 
resistance state is formed, Eq. (19), and the spontaneous 
dc field in the domains, Eq. (21). 

Our results are in overall agreement with the exper- 
imental findings. The observed T dependence of 
the photoresistivity at maxima compares well with the 
predicted T"^ behavior. Typical parameters uj^Tx ~ 
50 - 100 GHz, Tq ~ 10 ps yield cjrq/27r ~ 0.5 - 1 (over- 
lapping LLs), and the experimental data indeed closely 
resemble Fig. 2. For T ^ 1 K and ^ 100 K we find 
T;~^ ~ 10 mK, much less than ^ IK, as assumed in 
our theory. For the microwave power P ~ 1 mW and the 
sample area S ^ \ cm^, Eq. (20) yields the dimensionless 
power v'^^ ~ 0.005 — 0.1 (the smaller value corresponds 
to a;/27r = 100 GHz, the larger one to w/27r = 50 GHz), 
where we used rtr = 10 mK and vp — 2 - 10*" cm/s. These 

values of agree with characteristic values for sepa- 
rated LLs (Fig. 4) but are noticeably less than the predic- 
tion for overlapping LLs (Fig. 2). The discrepancy can be 
attributed (at least partly) to the fact that the value of 
Tq used in the above estimate, which was extracted from 
the Shubnikov-de Haas experiments, is in fact masked 
by inhomogeneous broadening and thus is shorter than 
the actual value. Indeed, Tq found from the (experimen- 
tal) damping of the oscillations in pph (which are local 
and thus not affected by inhomogeneous broadening) ac- 
cording to Eq. (15) is several times longer. With this 
value of Tq the threshold microwave power V*, Eq. (19), 
needed for the emergence of the zero-resistance states, 



corresponds to P less than 1 mW, in conformity with the 
experiments. ■^■^ Finally, for Vu; — V* ^ V*, at T ~ 1 K 
(when Ttr/Tin ^ 1), and LUc/27r = 50 GHz the estimated 
dc electric field in the domains, Eq. (21), is found to be 
^dc ~ 1 V/cm. This is in agreement with the experi- 
mental data of Rcf. 7 where the voltage drop between an 
internal and an external contact (separated by 200 /xm) 
generated by the radiation in the absence of the drive 
current was of the order of 5 mV for lUc/2tt ~ 20 GHz. 
Assuming that the size of the domain is of the order of 
the system size, this yields ~ 0.25 V/cm. One sees 
that this value indeed compares well with our theoretical 
estimate ~ 1 V/cm, especially taking into account 
the ojI dependence of E^^ following from Eq. (21). 

Recently, a number of publications appeared that ex- 
tended our theory (main results of which were pre- 
sented in Ref. 13) in a variety of contexts: propaga- 
tion of surface-acoustic waves, photoconductivity of 
laterally- modulated structures, ''^ photoconductivity for 
B above the cyclotron resonance, local compressibility 
of irradiated samples."^* 

We thank R.R. Du, K. von Klitzing, R.G. Mani, 
J.H. Smet, and M.A. Zudov for information about the 
experiments, and I.V. Gornyi for numerous stimulat- 
ing discussions. This work was supported by the SPP 
"Quanten-Hall-Systeme" of the DFG, by NSF grants 
DMR02-37296, DMR 02-13282 and EIA02-10376, by 
AFOSR grant F49620-01-1-0457, and by the RFBR. 



APPENDIX: DERIVATION OF THE BASIC 
EQUATIONS FROM THE QUANTUM 
BOLTZMANN EQUATION 

The purpose of this Appendix is to demonstrate how 
the semiclassical transport theory of Ref. 11 (see Sec. HI 
of this reference) is reduced to Eqs. (1),(3),(11) when the 
effect of electric fields on the impurity collision process 
can be neglected and only the distribution function is 
affected. The conditions under which the effect of the 
fields on the collision integral is weak are that the dc 
field is much smaller than Eg and the strength of the 
ac field satisfies P <C 1, where Eq and P are defined in 
Eqs. (5.5) and (6.2) of Ref. 11, respectively. In notation 
of our Eq. (10) these conditions read^® 



'Pu>, Qdc < Tin/Tq 



(Al) 



Under the conditions (Al) we neglect the effect of the 
external field on the density of states, which amounts to 
putting /ii = 1 in Eq. (3.42) of Ref. 11. The resulting 
DOS is time independent: 



1 -F2Re^A'ffj exp 



i2iTle 



(A2) 



where the coherence factor A = — exp(— Tr/wcTg), 
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and L[" is the Lagucrrc polynomial. The combination where 



|g;A'p has the meaning of the probability for an electron 
to complete the cyclotron orbit after / revolutions. Here 
we follow the notations of Ref. 11; in the main text of 
the present paper the factor A is denoted as —d. 

The kinetic equations (3.46) of Ref. 11 are written 
in the time representation for the distribution function 
f{t, t'; (fi, R), where (j) is the angle on the cyclotron orbit 
and R is the position of its guiding center. This function 
is related to the conventional distribution in the energy- 
time representation by the Wigner transform 



f{t,t 



, f de 



-ie{t-t') 



f 



t + t' 



(A3) 



The inverse Wigner transform of both sides of Eq. (3.46a) 
of Ref. 11 yields the canonical form of the Boltzmann 
equation for the distribution function 

{dt + uJcd^)f {t, e; if, R) = SUMe,t + Sti„{/}e,t. (A4) 

A time independent version Stin{/}e,i of the inelastic 
term is considered in Sec. VII and will not be discussed 
further in this Appendix. The impurity collision term 
Stim{/}£,t is obtained as the inverse Wigner transform 
of Eq. (3.46b) of Ref. 11 and will be written explicitly 
below. 

Under the condition (Al), the dissipative electric cur- 
rent at the point r is obtained from Eq. (A3) and 
Eq. (3.45b) of Ref. 11: 

jW(r,i)=2ep^ J gi(0) J ^J/(i,s;^,r,) 



+2Re^ A'^zexp ^^^^ f (t - ■^,s;V>,'rg 



(A5) 



where is the Fermi momentum, Vg = r^Rcei{(p)~C(t) 
is the guiding center coordinate and <^{t) describes the 
motion of an electron in the external field E(f): 



dtm 



eE{t) 



m 



(A6) 



i = {cos ip, sin (fi}, e is the antisymmetric tensor, e^y = 

^yx ~ 1* 

In view of Eq. (Al) one can neglect the effect of the 
electric fields on the interference processes in the collision 
integral (3.46b) of Ref. 11, which amounts to putting the 
form-factors /n = /12 = 1 in Eq. (3.49). Then Eq. (3.46b) 
in the time representation acquires a simple form 



1 °° ( 



Ttr 



(A7) 



Mi{t)-L{t,t') 



f; X'gi ( \M,{t') - L {t, t')] f (t, t' - ^) 



L{t, t') = {{ pF[C{t) - C(i') ]e - iRc^R } ■ i(<p) + id^ f 

and Mi{t) is the result of action of the operator 

L (t, t — 2-kI/lOc) on unity. Performing the Wigner trans- 
formation of Eq. (A7) we obtain 



Sti„,{/} 



1 / ' ' \ 2 
L(t+l-d^,t-l-de] fit, s; ip, R) + —Re V X^gi 



Mi(t + -d, 



X exp ( / ( i - — ,£;<p,R 



(A8) 



Equations (A4) and (A8) are valid for an arbitrary time 
dependent distribution function. We now notice from 
Eq. (A4) that the most divergent terms in the distribu- 
tion function are due to the contribution of the angular 
and time independent part of / to the collision integral. 
For WcTtr » 1 it is sufficient to consider / to be time, 
angular, and coordinate independent. Equation (A8) re- 
duces to 

St~{f}Z = -^^^f{e) + -Re V A'ffz 



N[^]-N{td,) 



(A9) 



and the bar stands for the time averaging. In the con- 
stant electric field N (St) = {5tYp\{dtCY /2, whereas for 
the microwave field C{t) = Re (Ca-e*"^*) and N {5t) = 
PfRs [CC(l-e*'^''*) ]/2- Substituting these expressions 
in Eq. (A9), taking C, from Eq. (A6), and using Eq. (A2), 
we arrive at the left-hand side of Eq. (9). 

The zero angular harmonic of the distribution func- 
tion does not contribute to the electric current (A5) di- 
rectly. The relevant angular dependent correction df{ip), 
with / d(pSf{(p) — 0, can be found perturbatively from 
Eqs. (A4) and (AS): 



l-t-2ReVA'5zexpf^') dj{e). (AlO) 
1^1 \ J _ 



Substituting the solution of Eq. (AlO) in Eq. (A5) and 
applying Eq. (A6) in the limit of the dc field, we obtain 
Eq. (1) with crdc(e) given by Eq. (3). 
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